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Abstract 

We report on a search for D^—D^ mixing made by a study of the 'wrong-sign' 
process ^ K^tt^ . The data come from 9.0 fb^^ of integrated luminosity 
of e~^e~ collisions at y/s ~ 10 GeV, produced by CESR and accumulated 
with the CLEO II. V detector. We measure the time-integrated rate of the 
'wrong-sign' process K~^7r~, relative to that of the Cabibbo-favored 

process W^K+tt', to be i?ws = (0.34 ± 0.07 ± 0.06)%. By a study of that 
rate as a function of the decay time of the D^, we distinguish the rate of direct, 
doubly-Cabibbo-suppressed decay ^ K^tt~ relative to ^ K~^tt~ , to be 
Rd = ± 0.08)%. The amplitudes that describe D^-W mixing, x' 

and y', are consistent with zero. The one-dimensional limits, at the 95% C.L., 
that we determine are (l/2)x'2 < 0.05%, and -5.9% < y' < 0.3%. All results 
are preliminary. 
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Neutral particles such as the K^, D^, B^, and mesons can evolve into their respective 
antiparticles, the K'^, D^, 5° and Measurements of the rates and mechanisms of 

K^—K^ and B^ — B'^ mixing have guided the form and content of the Standard Model, and 
permitted useful estimates of the masses of the charm and top quark masses, prior to direct 
observation of those quarks at the high energy frontier. 

Within the framework of the Standard Model, the rate of — £)° mixing is expected to 
be small, for two reasons. First, the decays of the are Cabibbo favored, while the mixing 
amplitudes are doubly-Cabibbo-suppressed, and second, the GIM cancellation is thought 
to cause substantial additional suppression. Many interesting extensions to the Standard 
Model predict enhancements to the rate of D^ — D^ mixing 0, much as the charm and top 
quarks enhance the rates of K^ — K^ and B^ — B^^ mixing. 

We report here on a study of the process, K^n' (those processes obtained by the 

application of charge conjugation on all particles, such as, in this case, D° K~n^, are 
implied throughout this report). We use the charge of the 'slow' pion, 7r+, from the decay 
Z}*+^Z}'^7r+ to deduce production of the D^, and then we seek the rare, 'wrong-sign' i^+vr" 
final state, in addition to the more frequent 'right-sign' final state, K^n^. 

The wrong-sign process, D'^ —*K~^7r~, can proceed either through direct, doubly-Cabibbo- 
suppressed decay (DCSD) via the Feynman diagram portrayed in Fig. |l|, or through mixing 
followed by the Cabibbo-favored decay (CFD), ^ ^ K'^n' . Both processes contribute 
to the 'wrong-sign' rate, 



To disentangle the two processes that could contribute to D'^ JC+vr", we study the 
distribution of wrong-sign final states as a function of the proper decay time, t, of the 
The mixing amplitude grows, relative to the decay amplitude, by a factor of — {1/2) {ix + y)t, 
as portrayed in Fig. ^. We refer to the proper decay time t in units of the mean lifetime, 
r£)0 = 415 ± 4fs [|]. Then, we refer to the mixing amplitude for — > D^, in appropriate 
units: those of one-half the mean decay rate, Tdo/2 = 1/[2tdo]. The mixing amplitude 
through virtual intermediate states is x, and that through real intermediate states by y [Q. 
When conservation of CP is assumed: 



X 



y 



AM _ Mi2 

Ar ^ ri2/2 

2T£)0 Vdo/2 



We use the convention that ?/ > 0(?/ < 0) corresponds to a shorter (longer) than average lifetime 
for eigenstates of CP with the same eigenvalue as the tt+tt" state populated by ^ti^tt~ . 

It is likely, although not inevitable, that within the Standard Model Exten- 

sions to the Standard Model contribute to x alone. 

The possible interference between direct decay and mixing would cause the number of 
wrong-sign decays, relative to the total number of right-sign decays, as a function of t, r„s(t) 
to be |igjTl[ 
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CcTT 



u u u 

c ^ d s c ^ s "d 

u u u 

B oc -e-^*^ sin^ Oc A (x cos^ 6>c 

= I^MP ~ tan"^ 6>c ~ 1/400 

FIG. 1. The Feynman diagrams for the decays Kir. The 'wrong-sign' decay, with am- 

phtude denoted B, is shown on the left, and the 'right-sign' decay, with amphtude denoted A, 
is shown on the right. The Cabibbo angle is 6c, and 5 is a possible relative strong phase be- 
tween the processes shown. The ratio of decay rates is denoted Rd and is expected to be of order 
tan4 6ic~ 1/400 = 0.25%. 
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rUt) = \[^-\{ix + y)t]e-'/^ 



= [Rd + ^/Roiycosd - xsm6)t + -{x^ + y'^)t'^]e * 
= [Rd + VRE y' t+^{x'^ + y'^)t^]e-' 

FIG. 2. The two processes that yield wrong-sign Kit. The proper decay time of the D^, in 
units of the mean hfetime, is denoted t. The strong decay of the D*~^ —>■ D^n^ specifies t = 0, 
and the charge of the slow pion 7rj" distinguishes a from a D^. The can undergo direct, 
doubly-Cabibbo-suppressed decay (DCSD), to K~^7r~ , schematically represented by the route to the 
left, with amplitude Se"*/^. The 1)° might also mix, with amplitude — (1/2) (ix + y)t, schematically 
represented by the route to the right, into a D^, and then undergo Cabibbo-favored decay (CFD) 
into K~^7r~ . The interference of these two amplitudes gives the relative decay rate to K~^7r~ as 
a function of t, r^s 

(t), shown. The quantity y' is defined to be ycos6 — xsinS, and x' is defined 
analogously, x' = x cos S + y sin 6. 
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rUt) = [Rd + \fR^y't + + y'^)f]e-\ (1) 

where Rd is the ratio of wrong-sign to right-sign rates of direct decay. The D^ — D^ mixing 
amphtudes x' and y' are defined by: 

y' = ?/cos(5 — xsin(5 
x' = X cos 5+?/ sin (5 ' 

where 5 is a possible strong phase between the wrong-sign and right-sign decay amphtudes. 
In a sense, use of a hadronic final state, such as Ktt, 'filters' the mixing amplitudes, much 
as a polarizer filters light. There are plausible arguments that y' > [0, and 6 < 13° 

We report here on the analysis of data accumulated between 1995 and 1999 from an 
integrated luminosity of 9.0 fb~^ of e'^e~ collisions with ^/s ~ 10 GeV at the Cornell Electron 
Storage Rings (CESR). The data were taken with CLEO II multipurpose detector |T2| , 
which includes two cylindrical drift chambers in a superconducting solenoid {B = 1.5 T) 
for measurement of the three-momentum p of charged particles, a cylindrical array of Csl 
crystals for measurement of the energies of photons and electrons, and a system of iron 
absorbers interleaved with proportional chambers for the identification of muons. These 
systems cover approximately 80% of the solid angle around the e~^e~ annihilation point. 
In 1995 a silicon vertex detector (SVX) was installed |T3|, that enables both precise re- 



construction of the proper lifetime of short-lived particles such as the D'^, as well as improved 
reconstruction of the angle 6 between the e~^e~ beams and the momentum vector of low- 
momentum charged particles, such as the slow charged pion from the decay D*~^ D^ir^. 
Also, in 1995, the gas in the larger CLEO II drift chamber was changed from argon-ethane 
to a helium-propane mixture, resulting in improved momentum and mass resolution, and 
improved particle identification by specific ionization (dE/dx). We refer to this revised 
configuration of the CLEO II detector as CLEO II. V. 

We reconstruct candidates for the decay sequences D*^ followed by either 

—y K~^'ir~ (wrong-sign, or WS) or —>■ K~n^ (right-sign, or RS). The sign of the slow 
charged pion, either 7r+ or vr", identifies ('tags') the charm state at production ('t = 0') as 
either or D^. The broad features of the reconstruction are similar to those employed in 



the recent CLEO measurement of the D meson lifetimes ||1^, but there are four principal 
differences. First, we accept candidates with total momentum, pd*, as low as 2.2 GeV; we 
show the approximate D* production cross section as a function of po* in Fig. ^. Second, 
we fully exploit the three-dimensional tracking capability of the SVX. We reject candidates 
where the daughter tracks form poor vertices in three-dimensional space, and we sharpen 
the reconstructed direction of the 7r+ momentum, thereby improving our resolution for re- 
construction of Q, the (small) energy released in the D*'^ '^t^^ decay. We denote by M 
the reconstructed mass of the two daughters of the and by Mk-ktt the reconstructed mass 
of all three particles in the candidate. We reconstruct Q = {Mk-k-k — M — mn)c^ , where m.„ is 
the mass of a charged pion. Third, we require candidates to be well-measured in Q, and in 
M. Fourth, we require candidates to pass two 'kinematic' requirements, designed to suppress 
backgrounds from D'' — > tt+tt", ^ K~ , ^ > 2 bodies, and from cross-feed between 
WS and RS. For the first kinematic requirement, for ^ K^ti~ candidates we evaluate the 
mass M under the three alternate hypotheses D"— ^tt^tt^, ^ K^K^ , and ^7^~^K~ . If 
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any one of the three masses so computed falls within 30 MeV/c^ (approximately 4.5 a) of the 
mass, the K^tt~ candidate is rejected. A conjugate requirement is made for the 

RS decays. The second kinematic requirement rejects one of the two configurations of very 
asymmetric decay, with a requirement that the decay angle 9* of the hypothesized kaon, 6*, 
in the D° rest frame, with respect to the D° boost, satisfy cos 6'* < 0.8. This requirement 
removes candidates with slow pions from the decay. 

D*- Production at 10.58 GeV 




p^^. (GeV/c) 

FIG. 3. The D*^ production cross section at ^ = 10.58 GeV @,[|]. The D*^ at higher 
momenta result from the continuum production: e~^e~ -^cc, followed by hadronization of charmed 
quarks to a D*^, and these D*^ are present for all of our data. The component of D*^ at lower 
momenta are present in the 2/3 of the data taken with ^/s at the mass of the T(4s), and these softer 
D*^ result from the decay B — > D*^X. We accept D*^ with momenta greater than 2.2GeV/c; 
most of these are from continuum production. 



The CLEO drift chamber system allows reconstruction of the specific ionization (dE/dx) 
deposited by the passage of charged particles, permitting an independent assessment of the 
identification of a charged particle as either a vr or a i^. We require that the dE/dx be within 
three standard deviations (cr) of the hypothesis used in the reconstruction of M and Q; this 
is a loose, consistency requirement. We tighten dE/dx criteria only for the evaluation of 
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systematic errors. 

A total of 16126 ± 127 candidates for the right-sign decay pass all requirements; their 
distribution in Q and M is shown in Fig. ^ 

The population of the wrong-sign candidates is shown on a scatter plot of Q versus M 
in Fig. ^ In Fig. ^ several prominent background features are evident. The combination 
of K^TT~ with random slow pions, vr"^, causes the vertical 'ridge,' with M equal to 

the mass of the M^o. The phenomenon is sometimes called 'dilution', and we will refer 
to it as 'random vr^ + D^/D^\ A diffuse background is also evident, due to the random 
combination of charged particles from e~^e~ flight quarks, and from e~^e~ cc, which we 
will refer to as combinatoric 'wrfs' and 'cc,' respectively. There is also a broad enhancement 
near the signal in Q, but at M < Mj^o, due to processes such as , followed by a 

very asymmetric decay p+ ^7r°7r+, resulting in a 7r° that is nearly at rest. Since a number 
of modes that we model can produce such behavior, including ^ p'^7r~ , D'^^p~n~^, and 
D'^ ^ K*^TT^ , we refer to this type of background as 'PV, for pseudoscalar-vector. 

To deduce the wrong-sign rate, R^^, we perform a 2-dimensional fit to the region of the 
Q versus M plane shown in Fig. |^. Each of the background shapes: 'random vr^ + D^/D^\ 
cc, uds, and PV is taken from Monte Carlo simulated data, which statistics corresponding 
to 90 fb~^ of integrated luminosity, but the normalization of each component is allowed to 
float in the fit; only the uds shows a significant difference with its expected contribution; 
the fit exceeds expectation by approximately a factor of two. The wrong-sign signal shape 
is taken directly from a seven a ellipse around the right-sign signal. The fitted background 
composition, superposed on the data, and projected on to Q and M are shown in Fig. |[ 
The event yields in the signal region from the fit are summarized in Table I. 

TABLE 1. Event yields in a signal region of 2.4 a centered on the nominal 
Q and M values, for the various categories of signal and background. The 
total number of candidates is 107. The first five rows are from the fit, in 
the two dimensions of Q and M, to the data shown in Fig. ^. Projections 
of the background components of the fit are shown in Fig. ^ The errors are 
statistical alone, and for the background components, are the errors on the 
mean yield in the signal region. The sixth line results from a fit to the data 
right-sign data in Fig. ^ 



Component ^ Events 

D^-^K+TT' (WS Signal) 54.8 ± 10.8 

random 7r± + D^W 24.3 ±1.8 

cc 12.3 ± 0.8 
uds 8.6 ± 0.4 

PV 7.0 ±0.4 



po^K+n- (RS Normalization) 16126 ± 126 

No acceptance corrections are need to directly compute, from Table I, i?ws = (0.34 ± 
0.07)%. The dominant systematic errors all stem from the potentially inaccurate modeling 
of the initial and acceptance-corrected shapes of the background contributions in the Q-M 
plane. We assess these systematic errors by substantial variation of the fit regions, dE/dx 
criteria, and kinematic criteria; the total systematic error we assess is 0.06%. 
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16,126 ± 127 Signal 
Gq = 190 ± 6 keV 
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FIG. 4. Signal for the right-sign process K~it~^ . The events in the top and bottom 

plots come from the same sample, and are all candidates for the decay D*^ T^t^^ followed by 

^ K~Tr^ . In the top plot, the horizontal axis is Q, the energy released in the decay of the D*^, 
for events where the reconstructed mass M of the Kir system falls within 15 MeV/c^ of the 
mass. In the bottom plot, the horizontal axis is M, for events where Q falls within 450 keV/c^ of 
its nominal value. The total number of right-sign signal events is 16126 it 127. The experimental 
resolution function in Q is non-gaussian, due to a substantial dependence of the width of the specific 
(gaussian) resolution on the particular decay configuration; the aq shown is the core value from a 
fit to a bifurcated Student's t distribution. The result for au comes from a fit to a double gaussian. 



10 



> 
o 



0.01 



0.02 F 







t .. :•.: •.- 



' ' ' 



1.8 



1.85 



1.9 



1.95 



M(GeV/c ) 



FIG. 5. Scatter plot of Q versus M for the candidates for D K^tt^ . The signal region is 
lightly shaded. The most prominent feature is the vertical 'ridge' at AI^o ; this results from the 
combination of the CFD ^ K^it~ with a random slow pion vr''", which combine to fake a £)*"*". 
The pervasive, smooth background consists of random combinations of charged tracks, referred to 
as cc and uds background in the text. The background density for M < M^o exceeds that for 
M>M£)0 simply due to phase space for the random combinations, and also due to true decays 
where a particle has been missed, resulting in a reduction in M. The so-called 'PV background, 
described in the text, causes an increase in density near the signal in Q, but at M < Mj^o. 
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FIG. 6. Signal for the wrong-sign process D^K^ir^. See also the caption for the right-sign 
analog, Fig. ^. The solid lines are the data, and the colored regions are the contributions from 
various background sources from the fit, as labeled, with shapes from simulation, and levels from a 
fit to the two-dimensional plane of Q and M. The fit to the signal peak is not shown. The results 
of the fit are summarized in Table I. 
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Our complete result for R^g is summarized in Table II. 



There are two directly comparable measurements of R„s- one is from CLEO II \\16 
Rws = (0.77 ±0.25 ±0.25)% which used a data set independent of that used here; the second 
is from Aleph |jT^, _R„s = (1-84 ± 0.59 ± 0.34)%; comparison of our result and these are 
marginally consistent with = 6.0 for 2 DoF, for a CL. of 5.0%. 

TABLE II. Result for R^s- For the branching ratio B{D^ ^ K^n^), we take 
the absolute branching ratio B{D^ ^ K~^-k~) = (3.85 ± 0.09)%, and the third 
error results from the uncertainty in this absolute branching ratio. 

Quantity Result 

R^g (0.34 ±0.07 ±0.06)% 

R^Jtan^Oc (1.28 ± 0.25 ± 0.21) 

B{D°^K+n-) (1.31 ±0.26 ±0.22 ±0.03) X 10"^ 

We have split our sample into candidates for —>■ K^tt^ and K^n^. There is 

no evidence for a CP- violating time-integrated asymmetry. From Table II, it is straightfor- 
ward to evaluate the la statistical error on the CP violating time-integrated asymmetry as 
v/T07/54.8 = 19%. 

Given the absence of a significant time- integrated CP asymmetry, we undertake a study 
of the decay time dependence of the wrong-sign rate in which CP conservation is assumed. 
Our fits use Equation |I], which describes the wrong-sign decay time dependence, include the 
term that is linear in t, unless specifically noted. The fit variable y' is allowed to vary over 
all real values, and we thereby account for the principal objection |18| to an earlier analysis 



of D^ — D^ mixing |]T9|. We believe that the fact that our acceptance in t extends all the 
way to zero lifetimes, while the acceptance for the fixed target experiments E691 and E791 
T9| , p!8[| dropped near t~l/2, prevents the loss in sensitivity that E791 |18| reported. 



We reconstruct t using only the vertical, or y, component of the flight distance of the D^. 
This reconstruction is guided by the physical dimensions of the CESR luminous region PD| 



and is portrayed in Fig. ^. We reconstruct the decay point, (x^, y^, z^), with a resolution 
that is typically 40 fim in each dimension. We measure the centroid of the luminous region, 
{xb,yb,Zb) with suitable hadronic events in blocks of data that typically contain integrated 
luminosities of several pb~^. Runs where yb is poorly measured are discarded, and the 
dominant error on y^ — yb comes from y^. We reconstruct t as: 

^ ^ M ^ yv-yb 
Py ctdo 

where py is the ^/-component of the total momentum of the K~^tt^ system. The error at is 
typically 1/2, although when the D° direction is near the horizontal, at can be large; we 
reject candidates with at > 3/2. 

In the sample of right-sign, K^-n^ decays, 14199 ± 120 signal events pass the 

requirements for the reconstruction of t, out of the initial sample of 16126 ± 126 used for the 
measurement of i?ws- The distribution in t of the right-sign events is shown in Fig. ^. We 
fit that distribution in the manner described in our measurement of mean charm lifetimes 
|1^, and find, in units of the recent world average @] mean life, 0.972 ± 0.014, where no 
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Reconstruction of t 



K+ 



+ 

I D*+ 




Py CTdo 



vertical 

(y) 



e+e beam 
annihilations 



CESR Plane (x) - 700 /im 



FIG. 7. The reconstruction of t. We measure the proper decay time of the D^, t, in multiples 
of the mean lifetime, or i = t/tj^o, where r is the proper decay time measured in seconds. 
As shown in this diagram, the region of e~^e~ annihilation is very small in the vertical, or y, 
where the region extends only 20 ;um at FWHM. In the bend plane of the CESR collider, or x, 
the FWHM extends 700 ^m, and along the beam axis or z, the FWHM is 30, 000 fj,m. The typical 
three-dimensional flight distance of the is 180 /xm, which is small compared to the x and z 
extents of the luminous region. We therefore reconstruct t using the displacement of the decay 
point, or vertex, in y from the beam centroid using the formula as shown. The symbol y^ refers 
to the K~^7r~ vertex, and yb refers to the beam centroid. The typical resolution in t is at ^ 1/2, 
and at becomes large when the direction of the flight is near to the horizontal. The error in 
extrapolation of the irf back to the luminous region is large, due to multiple scattering in the beam 
pipe, and indeed, reconstruction of the vr^ direction benefits greatly by constraint of the vr^ to 
the intersection of the beam and the extrapolated flight path. Use of other tracks in the event 
typically do not improve reconstruction of the estimate of the D*^ production point beyond that 
shown. 
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systematic error for this measurement is assessed. Since we have not assessed a systematic 
error, this result should not be construed as a new result on tdo . However, possible systematic 
errors due to the reconstruction and fitting technique are limited by this result from the 
right-sign data. There are far fewer events in the wrong-sign data, and so this class of 
systematic errors is negligible compared to the statistical error of the fit to the wrong-sign 
data. 

The fit to the right-sign data determines the resolution function that we use for the fit 
of the wrong-sign data. We also use the mean lifetime from the fit to the right-sign data 
as the central value for various charm backgrounds that are present in the wrong-sign data. 

Our resolution function is displayed, as well as the effect of folding that function with 
each of the three functional forms, e~*, te~*, and t^e~*, that enter in the time-dependent rate 
of wrong-sign decays from Equation |I], in Fig. |^. For the forms te~* and t^e~*, the breadth 
of the distribution is dominantly from the functional forms themselves; a narrowing of the 
resolution function would not greatly help distinguish between them. 

In the sample of wrong-sign, ^K^-k^ decays, 91 events pass the requirements for the 
reconstruction of t, out of the initial sample of 107 events in the signal region described in 
Table I, and used for the measurement of i?ws- The distribution in t of the 91 wrong-sign 
events is shown in Fig. |10[ We fit that distribution in the manner distinct from that used in 



our measurement of mean charm lifetimes 



We estimate the composition of the 91 events by applying the same requirement to the 
(renormalized) Monte Carlo sample used to fit for i?„s. A summary of the composition 
appears in Table III. 

We fit the distribution of wrong-sign candidates displayed in Fig. |10| with a superposi- 
tion of rws(t), from Equation |I|; an exponential with the lifetime for the 'vr^ -|- / D^^ 
background; a zero-lifetime component for the uds background; and second exponential to 
describe the cc and PV backgrounds. All distributions in t are folded with our resolution 
function. A Poisson likelihood, based on the extended maximum likelihood technique [pT| , 
is computed, for bins that are 1/10 of a lifetime, which is much smaller than at- 

The mean lifetime of the Monte Carlo simulated data, after renormalization for the fit to 
the Q—M plane, agrees very well with the data, outside of the signal region. This agreement 
gives us confidence in the use of the Monte Carlo simulated data to describe the background 
in the signal region. In the Monte Carlo simulated data, the decay time distribution of the cc 
and PV backgrounds is well described with a single exponential, folded with the resolution 
function. We use the mean D° lifetime to describe that exponential, but we explicitly vary 
that lifetime by ±15% to assess a systematic error. 

The information on the expected fractions of the various background components, as it 
appears in Table III, is incorporated as a a separate constraint added to the likelihood. The 
error on each fraction is nominally ((Jp + apY^"^, except for the signal component, where 
ap is omitted. The systematic errors, as from Table III, are included when the systematic 
errors on the fit parameters are evaluated. 
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FIG. 8. Distribution in t for right-sign D^^K-7r+. The histo gram shows the data, and those 
data are the same for the upper and lower plots; the only difference is the vertical scale, linear 
above, logarithmic below. The smooth curve is the fit, which results in the mean life shown of 
T£)0 = 403.5 ib4.1fs, where the error is statistical alone. We have not assessed the systematic error 
on the fit to r^o, so, it should not be construed as a new result on r^o. 
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FIG. 9. The functional forms of decay and mixing, as a function of the proper time t, relative 
to the mean lifetime. The top plot shows our resolution function, as determined from the fit 
in Fig. which can be approximately characterized by a gaussian with resolution at ~ 1/2. The 
subsequent three plots show this resolution function convolved with the three functional forms 
evident in Eq. |l|: e~*, te~*, and t^e~*, which describe, respectively, DCSD decay, interference 
between DCSD decay and mixing, and mixing. 
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TABLE III. Event yields in the Q — M signal region for wrong-sign K~^7t~ 
candidates that have had t successfully reconstructed. The distribution in t 
itself is shown in Fig. ^ The event yields in the second column are taken 
from the same fit to the data described in Table I, but with the effects of 
successful t reconstruction taken into account by the study of our sample 
Monte Carlo simulated data, which corresponds to an integrated luminosity of 
90 fb~^. The statistical errors on the number of events, in the second column, 
result from the fit to the whole Q — M plane; since the signal populates the 
signal region with high efficiency, the errors on the K~^n~ signal are 

the full Poisson errors on the mean yield of events. The statistical errors on 
the background components from the fit are relatively small compared to the 
Poisson error that corresponds to the mean event yield; a far greater number of 
background events populate the entire Q — M plane than populate the signal 
region alone. The third column gives the fraction of the wrong-sign (WS) 
sample constituted by the signal and each category of background; the fourth 
column gives the Poisson error ap, corresponding to any one random sample 
drawn from a distribution with the mean yield of events, on that fraction. The 
fifth and sixth columns give the error ap on each fraction from the fit to the 
whole Q—M plane, and then the error as from systematic effects, determined 
principally from varying fit regions in the Q—M plane. For combining errors, 
ap and ap are counted only once for the signal. The last column shows the 
mean life of the background components, with our estimated systematic error 
on that mean life, in units of the mean lifetime. For the signal ^ K^7r~ , 
the mean life is unknown (unk.) and can vary, according to Eqn. |I|, between 
0.586 and 3.414. For reference, the right-sign yield and measured mean life 
with error is given in the last row. 





# 


% 


% 


% 


% 




Component 


Events 


WS 


a-p 


o-p 


o"s 




WS Data 


91 


100 












49.7 ±9.8 


54.6 


10.8 


10.8 


8.7 


unk., .59-3.4 


Backgrounds: 














n^ + OyW 


19.6 ± 1.4 


21.5 


4.9 


1.5 


2.5 


.972 ±0.014 


cc 


10.4 ± 1.1 


11.4 


3.5 


1.2 


1.4 


1.00 ±0.28 


uds 


5.1 ±0.6 


5.6 


2.5 


0.7 


1.2 


.004 ±0.015 


PV 


6.2 ±0.7 


6.8 


2.7 


0.8 


0.8 


1.17 ±0.25 




14199 










.972 ±0.014 



Six parameters are able to vary in the fit: x', y' , ap where Rp = aj^, all from Equation |l|; 
and the fraction of each of the vr^ ± D^/D^, combined cc±PV, and uds backgrounds. 

The fit has been tested on 73 samples synthesized from the right-sign candidates, and 
background from the sidebands in the wrong-sign Q — M plane. Each of the synthesized 
samples contains 91 events, and they are composed to represent the wrong-sign sample 
described in Table III. The fits to the 73 synthetic samples show no bias in fitted values of 
x', y', and ap. They also verify the ability of the fit to estimate our statistical errors. 
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FIG. 10. Distribution in t for wrong-sign ^K^t:^ . The same histogram of data is shown in 
both the upper and lower plots; the differences are the fit results shown and overlayed. The solid 
black lines are the totals of all components of the fits. The values of x' and y', from Eqn. |l|, are fixed 
to zero in the fit displayed in the upper plot. The best-fit shapes from total and uds backgrounds 
are shown as smooth curves. The variables x' and y' float freely in the lower fit, and the best-fit 
values are shown. The likelihood decrement from the upper to lower fit gives V— 2Aln£ = 1.8 o". 
The best-fit shapes of the two mixing terms are shown, although only that from te"* is visible. 
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In the initial fit to the actual wrong-sign data, x' and y' are constrained to be zero, and 
this fit is shown in the upper part of Fig. [1^; the confidence level of that fit is 84%, indicating 
a good fit. 

The mixing amplitudes x' and y' are then allowed to freely vary, and the best fit values 
are shown in both the lower portion of Fig. and in Table IV. 

TABLE IV. Resuhs of the fit to the distribution of D^^K+ti' in t. Both the 



distribution, and the fit, are shown in the lower portion of Fig. ITO. 



Parameter 


Best Fit 


95% C.L. 


Rd 

y' 


(0.52t°-l^ ± 0.08)% 
(-2.7i||±0.2)% 


-5.9% <y' < 0.3% 


x' 

(l/2)x'2 


(0± 1.6 ±0.2)% 


|x| < 3.2% 
< 0.05% 



The fit improves slightly, by an amount corresponding to \^^2A\nC = 1.8a, including 
systematic effects, when mixing is allowed. Our interpretation of this change is that it 
represents a statistical fluctuation. 

Therefore, our principal results concerning mixing are the one-dimensional intervals, 
which correspond to a 95% confidence level, that are given in Table IV. 

Additionally, we evaluate a contour in the two-dimensional plane of y' versus x', which 
at 95% confidence level, contains the true value of x' and y'. To do so, we determine the 
contour around our best fit values where the — ln£ increases by 3.0 units. The interior of 
the contour is shown, as the small red region, near the origin of Fig. 0. On the axes of x' 
and y', this contour falls slightly outside the one-dimensional intervals listed in Table IV, as 
expected. 

We have evaluated the allowed regions of other experiments , at 95% C.L., 



and shown those regions in Fig. 0. In combining the E691 and E791 studies of hadronic 
final states, we make the most optimistic assumption, that leads to the smallest allowed 
region, concerning treatment of the term linear in t in Equation |l|. We do not necessarily 
endorse that interpretation, and we note that the E791 results utilizing D'^ K^i~ui suffer 
no similar uncertainty in interpretation. 

Finally, if we assume that S is small, which is plausible [^^, then x' ~ x, and we can 
indicate the impact of our work in limiting predictions of D^ — D^ mixing from extensions 



to the Standard Model. We plot our one-dimensional allowed region in x' in Fig. |12|. Eigh- 
teen of the predictions tabulated in a paper contributed to this symposium have some 
inconsistency with our limit. Among those predictions, some authors have made common 
assumptions, however. 

In conclusion, we conducted a study of the wrong-sign process, ^ K^n^ , and con- 
clusively established its rate, relative to the right-sign process, D° K^tt~, as R^s = 
(0.34 ± 0.07 ± 0.06)%. By a study of that rate as a function of the decay time of the D^, 
we distinguish the rate of direct, doubly-Cabibbo-suppressed decay D^-^K'^t:~ relative to 
W-^K+n-, to be Rd = (0.50l|j:l^ ± 0.08)%. The amplitudes that describe D'^-W mixing, 
x' and y', are consistent with zero. The one-dimensional limits, at the 95% C.L., that we 
determine are (l/2)x'^ < 0.05%, and —5.9% < y' < 0.3%. The limit on x\ combined with 
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FIG. 11. Limits in the y' v. x' plane. Our experiment limits, at 95% C.L., the true values of x' 
and y' to occupy the dark (red) region near the origin. Also shown are the similar zones from other 
recent experiments. We assume (5 = to place the recent work of E791 that utilized ^ K~; 
a non-zero 5 would rotate the E791 confidence region clockwise about the origin by an angle of 6. 
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FIG. 12. Predictions impacted by this work. The horizontal axis is the Reference Index used 
in Q, that corresponds to the reference tabulated there. The horizontal region (red) nearest the 
bottom is the allowed region, at 95% CL, consistent with this work, in the plausible approximation 
that x' = x. There are 18 predictions that have some inconsistency with this work. 
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the plausible assumption that the relative strong phase 5 = 0, has some inconsistency with 
a variety of extensions to the Standard Model. Limits in the two-dimensional plane of y' 



versus x' are given in Fig. ITT 



All results described here are preliminary. 
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